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Abstract 
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io Using density functional theory plus Hubbard U calculations, we show that the ground state of 

§ ■ ii (Mg,Fe)(Si,Fe)03 perovskite, the major mineral phase in the Earth's lower mantle, has high-spin 

wo: 

^ . 12 ferric iron {S = 5/2) at both dodecahedral (A) and octahedral (B) sites. With increasing pressure, 

o ; 

13 the B-site iron undergoes a spin-state crossover to the low-spin state (S = 1/2) between 40 and 

14 70 GPa, while the A-site iron remains in the high-spin state. This B-site spin-state crossover is 



15 accompanied by a noticeable volume reduction and an increase in quadrupole splitting, consistent 

16 with recent X-ray diffraction and Mossbauer spectroscopy measurements. The anomalous volume 

17 reduction leads to a significant softening in the bulk modulus during the crossover, suggesting a 



00 ■ 

ib possible source of seismic-velocity anomalies in the lower mantle. 



19 PACS numbers: 91.60.Pn, 76.80. +y, 91.60.Gf, 91.60.Fe 
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20 The total electron spin (S) of a transition-metal ion in a crystalline solid can change with 

21 many factors, such as pressure, strain, or temperature, to name a few. This phenomenon, 

22 known as spin-state crossover, is of great importance in spintronics, as it allows artificial 

23 control of magnetic properties of materials, including coordination complexes with potential 

24 for molecular switches Not as widely known, spin-state crossover also plays a crucial role 

25 in geophysics. A well studied example is ferropericlase, (Mg,Fe)0, the second most abundant 

26 mineral (~ 20 vol%) in the the largest single region (~ 55 vol%) of the Earth's interior - 

27 the lower mantle. With increasing pressure, ferrous iron (Fe 2+ ) in this mineral undergoes a 

28 crossover from high-spin (HS) state, S — 2, to low-spin (LS) state, S — 0, in the pressure 

29 range of 40-55 GPa (flfl. The intermediate-spin (IS) state, S = 1, is not observed in 

30 this mineral. The HS-LS crossover in ferropericlase directly affects the structural, elastic, 



31 optical, and conducting properties of this mineral |6Nl lj and thus affects mantle properties. 

32 0,Q,Q|- 

33 In contrast, the spin-state crossover in iron-bearing magnesium silicate (MgSiOa) per- 

34 ovskite (Pv), the most abundant mineral (~ 75 vol%) in the lower mantle, has been a source 

35 of controversy for two main reasons. One is the coexisting ferrous and ferric iron (Fe 3+ ) in 

36 this mineral with an imprecisely estimated population ratio; the other is the lack of definitive 

37 tools to directly probe iron spin state at high pressures. Two techniques, X-ray emission 

38 spectroscopy (XES) and Mossbauer spectroscopy, have been widely used, but their inter- 

39 pretation can be ambiguous. The very similar XES spectra and M 6 ssbaue r spectra 



[16l-ll9| have been interpreted in terms of HS-IS and HS-LS crossover in (Mg,Fe)Si03 Pv. 
Plenty of calculations on (Mg,Fe)Si03 Pv have been conducted 0-0, but consistency with 



experiments was not achieved until very recently 2j3|,|27j]. Now the spin state in (Mg,Fe)Si03 
Pv is better understood: the observed increase of iron nuclear quadrupole splitting (QS) in 
Mossbauer spectra results from neither HS-IS nor HS-LS crossover, but from the change in 
the 3c? orbital occupancy of the HS iron 27J. As to ferric iron in Pv, possibly more abundant 
than ferrous iron ( Fe 3+ / J^Fe might be as high as 2/3) 28|, |29], its spin-state crossover has 
remained unclear, as described below. 



48 Previous experiments investigating t 



re iron spin state in aluminum-free MgSiC>3 Pv were 



49 focused mostly on ferrous iron 16|, |l8|. Nevertheless, it was still observed that the low 



concentration of ferric iron in the sample exhibited an increase in QS with pressure, which 



51 suggests a crossover from HS (S = 5/2) to LS (S = 1/2) state in the pressure range of 

52 30-70 GPa. In contrast, in Al-bearing samples, where ferric iron occupies the dodecahedral 

53 (A) site, the QS remains unchanged up to 100 GPa, which suggests the A-site iron remains 

54 in the HS state [l7|. These results indicate that the ferric iron at the octahedral (B) site 

55 undergoes a spin-state crossover. Such a mechanism was recently confirmed by experiments 

56 using (Mgi_ a; Fe a ;)(Sii_a;Fe a ;)03 Pv (x = 0.1) samples: about half of the HS iron changes to 

57 LS state in the 45-60 GPa range while the other half remain in the HS state all the way to 

58 150 GPa [30]]. So far, the computational studies on (Mgi_ x Fe x .)(Sii_ x Fe x )03 Pv have found 

59 a ground state with HS iron at the A-site and LS iron at the B-site (A-HS; B-LS) and an 

60 A-site HS-LS crossover that leads both A- and B-site iron to a final LS state (A:LS; B-LS) 



6i at high pressures 
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22J. These predictions are inconsistent with experiments in two ways: 

62 (1) the predicted transition pressure is too high; (2) the predicted HS iron concentration is 

63 too low. 



To compare with recent experiments 30(, we stabilize (Mgi_ x Fe x )(Sii_ x Fe x )03 Pv with 
x = 0.125 in all possible spin states using a 40-atom supercell shown in Fig. [TJ We also 
calculate the iron nuclear electric field gradient (EFG) associated with each state, as the nu- 
clear hyperfine interaction has proven to be a unique fingerprint to identify the spin states of 
transition-metal ions 
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3l| . The atomic structures were fully optimized with damped vari 



69 able cell shape molecular dynamics 



32] implemented in the quantum espresso code |33J], 



where the plane- wave pseudopotential method is adopted 34j. These states were also inde- 
)endently confirmed via the augmented plane-wave plus local orbitals (APW+lo) method 



35] implemented in the WIEN2k code 



with which the EFGs were calculated. The EFGs 



were converted to QSs with 57 Fe nuclear quadrupole moment Q = 0.16 37J and 0.18 barn for 
the possible uncertainty. To treat (Mg 1 _ x Fe x .)(Sii_ x Fe x .)03 Pv, the density functional the- 
ory plus Hubbard U (DFT+£/) method is necessary, as standard DFT exchange-correlation 
functionals, the local density approximation (LDA) and generalized gradient approximation 
(GGA), sometimes lead to unwanted metallic states (especially at high pressures), in which 
the iron spin states are not well defined. Since the Hubbard U of A- and B-site iron in each 
spin state is unknown, we have to stabilize the desired spin state with a trial U and then 
extract the self-consistent U, referred to as U sc , using the linear response approach [3^] in 
a recently developed iterative procedure. This procedure is equivalent to, but more efficient 



82 than the one published earlier |39l|. and has been successfully implemented 40J. More details 

83 are described in the EPAPS |41 |. 



84 Within DFT+U, several combinations of iron spin states can be stabilized. The A-site 
as ferric iron can be stabilized in HS, IS, and LS states. The B-site ferric iron can be stabilized 
86 not only in LS state, but also in HS state that has not found in previous calculations 
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22]. 



87 The spin moments of the A- and B-site iron can be either parallel or anti-parallel. The U sc 
ss of ferric iron in Pv, listed in Table I, mainly depends on the iron spin state, slightly depends 
89 on the occupied site, and barely depends on pressure and alignment of spin moments. 



90 The relative enthalpy {AH) of each stabilizec 



state is shown in Fig. [21 where the pre- 



viously perceived ground state (A-HS; B-LS) |2l|, |22| is used as a reference. Remarkably, 
the actual ground state of (Mg,Fe)(Si,Fe)03 Pv has HS iron on both sites (A-HS; B-HS), 
regardless of the choice of exchange-correlation functional (LDA or GGA) and Hubbard U 
(U sc or 4 eV). These choices do not affect the spin-state crossover either: an HS-LS crossover 
only occurs in the B-site iron, while the A-site iron remains HS. As expected, the predicted 
transition pressure (Pt) depends on the exchange-correlation functional and Hubbard U : 



with LDA+U SC , P T = 41 GPa; with GGA+U SC , P T = 70 GPa; with GGA+U 



[U = 4 eV), 



Pt = 29 GPa. (Coordination complexes also show similar dependence [42|, |43|.) Notably, 
the alignment of iron spins (parallel or anti-parallel), barely affects Pt, as shown in Fig.[2]^c). 
The Pt predicted by LDA+U SC and GGA+U SC best agree with the Pt observed in Mossbauer 
spectra, 50-60 GPa [30|. The LDA+C/ SC electronic density of states (DOS) of the two relevant 
states (A-HS; B-HS and A-HS; B-LS) can be found in EPAPS Q- 

The calculated QSs of ferric iron (A- and B-site) and ferrous iron (A-site) {2^ in various 



104 spin states, along with the measured QSs 



16 
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30], are shown in Fig. |3j Clearly, our 



105 calculations on ferrous and ferric iron in Pv are consistent with Mossbauer spectra. The 

106 HS-LS crossover in the B-site ferric iron also helps to explain the decrease in the XES 

107 satellite peak (K/3') intensity [l4 , - Interestingly, the QS of ferrous and ferric iron exhibit 

108 exactly the opposite trends with respect to the spin moment. This can be understood via 

109 their orbital occupancies. The LS ferrous iron, although occupying the A site, is effectively 
no located near the center of a Fe-0 octahedron, as it is vertically displaced from the mirror 



plane 



241 ] . Its six 3d electrons doubly occupy the three orbitals with t 2g character and form 



112 a charge density with cubic-like shape |24j, which barely contributes to the EFG and leads 

113 to a very small QS. The HS ferric iron also has a small EFG (and thus QS), irrespective of 
H4 A or B site. This is because its five 3d electrons (all spin-up) occupy all 3d orbitals, forming 
us an almost spherically shaped electron charge distribution that leads to a small EFG (and 
lie thus QS). Similarly, the spin-up electrons in HS ferrous and LS ferric iron barely contribute 
in to EFG, as their charge distributions are nearly spherical and cubic, respectively. It is their 
us spin-down electrons that contribute to the EFGs and lead to larger QSs. This is why the 
n9 spin moments of ferrous and ferric iron appear to affect the QSs in an opposite manner. 



The LDA+U SC compression curves and bulk modulus (K = —VdP/dV) of (Mg 1 _ x Fe rr )(Sii_ 2; Fe :r )0 
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121 Pv (x = 0.125) along with the experimental data (x = 0.1) 30j are shown in Fig. HJ At 

122 low pressures (< 45 GPa), the experimental data falls on the calculated compression curve 

123 corresponding to the (A-HS; B-HS) state. Starting from ~45 GPa, the data points deviate 

124 from the (A-HS; B-HS) curve and then join the (A-HS; B-LS) curve at ~60 GPa. Starting 

125 from ~100 GPa, the data deviates from the curve again. This, however, is very likely to 

126 result from the questionable accuracy of the Au pressure scale used in the experiment, as 



127 already discussed in the case of (Mg,Fe)Si03 Pv 44|. Notice that the observed volume 

128 reduction further confirms the B-site HS-LS crossover, as the previously perceived A-site 

129 HS-LS crossover barely leads to a volume reduction, evident from the compression curves 

130 (A-HS; B-LS and A-LS; B-LS) shown in Fig. S(a). The B-site spin-state crossover and the 

131 observed volume reduction in the 45-60 GPa range can be qualitatively understood via the 

132 Fe 3+ electronic configurations and Fe-0 distances at A and B sites. With all 3d orbitals 

133 occupied, HS iron has spherically-shaped electron charge density and the largest radius 

134 compared with other spin states, favoring longer Fe-0 distances. Residing in the large do- 

135 decahedral cage, the A-site iron can easily maintain longer Fe-0 distances and thus remain 

136 in HS state. In contrast, the Fe-0 octahedron has smaller size and shorter Fe-0 distances. 

137 With increasing pressure, the internal octahedron bond lengths can be shortened enough to 

138 induce the HS-LS crossover. Since the 3d electrons of the B-site LS iron only occupy the 

139 t2g-hke orbitals pointing away from oxygen, the associated Fe-0 distances are significantly 

140 shorter than those of the HS iron at the same pressure. Therefore, the spin change of the 

141 B-site iron is accompanied by a noticeable octahedral (and thus unit-cell) volume reduction. 

142 Such volume reduction leads to anomalous softening in bulk modulus, as described below. 



143 At finite temperatures, the spin-state crossover passes through a mixed-spin (MS) state 

144 (namely, HS and LS coexist) within a finite pressure range that increases with temperature. 

145 During the crossover, the thermodynamic properties of the MS state exhibit anomalous 
we behavior that may affect mantle properties. One example is the softening in bulk modules 
147 and its effect on the compressional wave velocity, as already seen in ferropericlase 7J, |l(J, 111 ]. 



us To estimate such anomaly in (Mg,Fe)(Si,Fe)03 Pv, we employ a thermodynamic model 

149 similar to that used in Ref. (lOj. Here, we do not include vibrational free energy, as it 

150 barely affects the magnitude of the anomaly, slightly increases the transition pressure, and 



151 uniformly decreases the bulk modulus, as shown in the case of ferropericlase |10l.lll|. Indeed, 

152 the calculated V(P) curve of (Mg,Fe)(Si,Fe)03 Pv in the MS state (using LDA+U SC ) at room 

153 temperature (300 K), shown as the dashed line in Fig. HJa), exhibits a volume reduction 

154 (~1.2%) around the predicted Pt, 41 GPa. This reduction leads to a significant softening 

155 in bulk modulus, as shown in Fig. |4^b). The softening is still prominent at 2000 K, the 

156 temperature near the top of the lower mantle (~660 km deep). Given the abundance of 

157 iron-bearing Pv and the possibly high population of ferric iron, this softening may have a 
i5B noticeable impact on the mantle properties, including possible anomalies in the seismic wave 
159 velocities. 

wo In summary, with a series of DFT+U calculations, we have shown that the actual ground 

lei state of (Mg,Fe)(Si,Fe)0"3 perovskite has high-spin ferric iron on both A and B sites. It 

162 is the B-site ferric iron that undergoes a crossover from high-spin to low-spin state with 

163 increasing pressure, while the A-site iron remains in the high-spin state. The calculated 

164 quadrupole splittings and the compression curves are consistent with experiments. The 

165 volume reduction accompanying the B-site HS-LS crossover leads to a significant softening 
lee in bulk modulus, which suggests a possible source of seismic- velocity anomalies in the lower 
167 mantle. This work, one more time, demonstrates that the nuclear hyperfine interaction, 
lea combined with first-principles calculations, can be a useful tool to identify the spin states of 

169 transition-metal ions in solids under high pressures. 
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220 Table I. U sc , the self-consistent Hubbard U (in eV), of ferric iron on the A and B site in 

221 each spin state. 
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223 Figure Captions 

224 Fig. 1. (Color online) Atomic structure of (Mg .875Fe .i25)(Sio.875F e o.i25)03 Pv, configured 

225 with the shortest iron-iron distance, viewing along the [001] direction. Large (orange) and 

226 small (green) spheres represent Fe and Mg sites, respectively. Si-0 and Fe-0 octahedra are 

227 shown in opaque (blue) and translucent (orange) colors. 

22B Fig. 2. (Color online) Relative enthalpies of (Mg .875Feo.i25)(Sio.875Feo.i25)03 Pv in differ- 

229 ent spin states obtained using different functionals and Hubbard U. The reference state has 

230 HS iron in the A site and LS iron in the B site (A-HS; B-LS). Predicted transition pressures 

231 by LDA+U SC (a), GGA+U SC (b), and GGA+U with U = 4 eV (c) are 41 and 70, and 29 

232 GPa, respectively. Dashed lines in (c) correspond to anti-parallel spins at A- and B-sites. 



233 Fig. 3. (Color online) Calculated QSs of (a) ferrous iron [27J and (b) ferric iron in MgSi0 3 



234 Pv. Letter A and B in (b] 

235 effect of pressure on QSs 



refer 



16 



;o iron-occupying site. Arrows in (c) indicate the measured 



18 



3Q| 



p 



Fig. 4. (Color online) Compression curves (a) and bulk modulus (b) of (Mg 1 _ 2 .Fe x )(Sii_ x Fe a; )03 
v computed with LDA+U SC (x = 0.125) and room-temperature measurements (x = 0.1) 



23B 30[. Both the measured and calculated compression curves exhibit a clear reduction accom- 

239 panying with the B-site HS-LS crossover, which leads to a softening in bulk modulus shown 

240 in (b). 
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